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I Introduction 



Completely integrable systems in classical mechanics are well known to be of great interest for both 
theory and applications, see e.g. [H [2] and references therein. Indeed, the possibility of analytical 
description of the corresponding dynamics enables us to uncover important physical properties of 
the systems in question. The prime examples of this are the Kepler laws and numerous physical 
models based on (superposition of) harmonic oscillators. 

Interestingly, a number of physically relevant exactly solvable models (e.g., the Coulomb prob- 
lem for the hydrogen atom and the multidimensional harmonic oscillator, to name just a few) 
are maximally superintegrable rather than just completely integrable. In general, a Hamiltonian 
dynamical system on a 2n-dimensional phase space is maximally superintegrable, if it possesses 
the maximal possible number, 2n — 1, of functionally independent, globally defined integrals of 
motion (contrast this with n commuting integrals of motion for completely integrable systems); 
see e.g. PJ H] and references therein for further details. In this case we shall also say that the 
Hamiltonian of the system in question is maximally superintegrable. 

Quantizing a generic completely integrable classical system and solving the resulting Schro- 
dinger equation may often represent a nontrivial problem. However, solving the Schrodinger 
equation for maximally superintegrable systems is often easier (sometimes to the extent of reduc- 
ing the determination of energy spectrum to a purely algebraic problem, as is e.g. the case for 
the nonrelativistic hydrogen atom and the multidimensional harmonic oscillator) because of the 
presence of additional integrals of motion [3]. 

These facts have lead to a considerable interest in superintegrable systems in general, and 
maximally superintegrable systems with natural Hamiltonians on two- or three-dimensional con- 
figuration space are now quite well understood, see e.g. [6]-[T6] and the survey [Ij. However, much 
less is known about superintegrable systems in higher dimensions, although such systems often 
can be interpreted as multiparticle systems consisting of several one-, two- or three-dimensional 
particles with nontrivial interactions and therefore also can be of interest in physics. 

Moreover, quantum systems on an 2n-dimensional phase space that are exactly solvable for ar- 
bitrary n are relatively scarce, with a few beautiful exceptions like the multidimensional Coulomb 
problem and its extensions [T7] and the Calogero-Moser-Sutherland and related models, see e.g. 
[T8] and references therein. Therefore, any new examples would be of considerable interest. Higher- 
dimensional superintegrable systems are prime candidates to yield such examples: for instance, 
the above examples are superintegrable, see e.g. [H [T71 [19] and references therein. In the present 
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work we show that for two large classes of maximally superintegrable Hamiltonians on an 2n- 
dimensional phase space with arbitrary n the exact solutions for the Hamilton- Jacobi and the 
Schrodinger equation are readily available, see Theorems [3] and H] below for details. 

More specifically, we show that solving the Hamilton-Jacobi and the Schrodinger equations 
for the systems in question on the n-dimensional configuration space with arbitrary n amounts 
(provided n is sufficiently large) to solving fc-dimensional reduced equations, where k is fixed and 
independent of n, as presented in Theorems [3] and |H Moreover, the said reduced equations can 
be readily solved, as illustrated by the example of an anharmonic oscillator in Section V. 

In order to establish these results, we reveal a hidden symmetry of the systems under study. 
This symmetry manifests itself upon passing to the fiat coordinates of the metric tensors associ- 
ated with the kinetic-energy parts of the Hamiltonians under study, as described in Theorem [2] 
below. Namely, some of these fiat coordinates are cyclic coordinates for the Hamiltonians under 
study, so the corresponding momenta are integrals of motion, and the separation of variables is a 
reduced one in the sense of [H]. The dynamics in the non-cyclic coordinates is described by the 
reduced equations mentioned above. 

Note that the metric tensors in question have the so-called maximally balanced signature (the 
numbers of plus and minus signs differ at most by one), and the above fiat coordinates are the 
light-cone ones rather than the orthogonal ones. This leads to an interesting phenomenon: the 
aforementioned reduced equations are first-order rather than second-order PDEs, which makes 
them easier to solve. 

II Preliminaries 

For fixed integer m and k consider the separable Benenti Hamiltonians [22l [23] on the phase 
space T*Q, the cotangent space of an n-dimensional Riemannian manifold Q endowed with the 
contravariant metric tensor Gm'- 



Here A = (A^, . . . , A")"^ are coordinates on Q and fi = (/ii, . . . ,/Xn)^ are the corresponding mo- 
menta. The contravariant metric tensors Gm have the form [24\ 



where Aj = ]^ (A* — A-'), and L = diag(A^, . . . , A") is a (1, l)-tensor on Q called a special conformal 
Killing tensor [25J. The Killing tensors Kr from ([T]) are diagonal in the A-coordinates and can be 



r = 1, . . . , n. 



(1) 
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constructed as follows 



r-l 



\ r = 2,...,n, 



(2) 



fc=0 



where I is an n x n unit matrix, and cr^ = cr,fc(A) are symmetric polynomials in the variables 
A^, . . . , A" (ctq = 1,0"! = X]r=i -^N ■ ■ ■ ) = A^A^ ■ ■ • A"). They are related to the coefficients of the 
characteristic polynomial of the tensor L as follows: 



det(P-L) = ^(-l)V,^-^ 

1=0 

Consider the (g, p)-coordinates defined as follows [2U] : 

n 

= (-1)V,(A), = - 5^(A'=)"- WAfc, ^ = 1, 

fc=i 



Notice that are nothing but the coefficients of the characteristic polynomial of L 
(g, p)-coordinates we have 

n-l 



i+l 
j ' 



whence for m = 0, . . . , n we find [26] 



s:r+s 
"n— m+1 



n— m— 1 

E 



E ^-''^n-m+f+i, r,s = l,...,n-m, 



j=n— m+l 

otherwise. 



n-m+j+i' r,s — n m + 1, 



,n, 



(3) 



(4) 
In the 



Here and below 6j stands for the Kronecker delta. 

The geodesic Hamiltonians Em,r = ^fJ'^KrGmfi' are polynomial in the {p, g)-coordinates for 
m = 0, 1, 2, . . . [20]. In particular, for m = 0, . . . , n we have 

n—m n—m—l n~m m k 

Em,l 2 PjPn—m—j+1 ~l~ 2 5-^ ^ ^—j PjPn—m+k—j-\-l 2 5-^ ^ 5-^ Pn—m+jPn—m+k—j+l' 

j=l k=l j=k+l k=l j=l 

The basic separable potentials are given by the recursion relations [20] 



y(°) = 0, 



K+\ + Vr^'>Vr, k=l,2,..., w 



(1) 



A; = 1, 2, ... , K 



(-1) 



(5) 



where we tacitly assume that g° = 1 and that g* = for i > n. 



For any fixed m and k the Hamiltonians Hr"^''^^ are in involution with respect to the canonical 
Poisson bracket on T*Q 

These Hamiltonians are automatically separable in the (A, /i)-coordinates because they satisfy the 
Stackel separation relations by construction p4j. Note that the transition from the (A,/i)- to the 
{p, g)-coordinates is a canonical transformation. 



Ill Superintegrability and flat coordinates 

For any natural number n > 2 let 

ki 

Hm,r — -^mjr ~l~ ^ ^ ^k^r \ ^ ^i ■ ■ ■ i '^i 

k=ko 

where Ck are arbitrary constants and m, ko and ki are integers. It is readily seen that these 
Hamiltonians are in involution for any fixed m, ko and ki (see Theorem 1 of [27j for details): 

{HjYi^r: HfYi g^ 0, r, s 1, . . . , ?T,. 

We have the following straightforward generalization of Theorem 1 from [20] (see also p7|): 

Theorem 1 Given a natural n > 2 and an integer m G {0, ... ,72 — 1}, let k^ = —m and ki = 
n — 1 — m. Then Pn-m is an integral of motion for Hm,i, i-e., {ifm,i,Pn-m} = 0. 

Moreover, Fm,s = {Hm,s,Pn-m} are additional integrals of motion for Hm,i: {Hm,i, Fm,s} = 0, 
s = 2, . . . ,n, and the {2n — 1) integrals of motion for Hm,i (Hm,r, r = 1, . . . ,n, and Fm,s, s = 
2, . . . ,n) are functionally independent, so Hm,i is maximally superintegrable. 

Under the assumptions of Theorem [1] the Hamiltonian Hm,i involves n — 1 parameters Cj (note 
that V}^^ = 0). 

Proposition 1 ([26j) Given a natural n > 2 and an integer m E {0, . . . ,n}, the metric Cm in 
the coordinates r* defined by the formulas 

i-l 

gi — _j_ i ^ ^j^i-j ^ i = 1, . . . ,n — m, 

n . . (6) 

= '^'^'^^ 2 = n — m + 1, . . . , n, 

j=i 

takes the form 
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The transition from {q,p)- to (r, n) -coordinates, where tt^ = ^ TTT^'^' ^ = 1; • • • ? ^■5 canonical 

i=l OT 

transformation, and we have 

^ /n—m n 

Em,l = = 2 I 7rj7r„_m+i_j + ^ 7rj7r2„_m+l- 

\ j=l j=n—m+l 

The tensor L in the coordinates r* takes the form: 

for m < n: L) = 5]-^^! - 5'^_J - - ^r''-^-m+i+n[(j+m-i)/n]§'^_^ 

for m = n: L] = S]'^^ + IrV""-'"'"^. 

Here and below [k] denotes the largest integer less than or equal to k. 

The canonical coordinates (r, vr) are not orthogonal, but the metric tensor Gm is constant in 
these coordinates. Bringing Gm into the canonical form, with +1 and —1 at the diagonal and 
zeros off the diagonal, is possible [2B] but we shall not need this here. 

Recall [20] that for k = 1, . . . ,n — 1 the potentials v}''^ are independent of with j = 

t—k) 

+ l,...,n. Likewise, for k = 1, . . . ,n — 1 the potentials are independent of with 

j = 1, . . . ,n — k. On the other hand, the change of variables (El) is partially triangular: with i = 
1, . . . , n — m depend only on r^, . . . , r* while with i = n — m + 1, . . . ,n depend only on r*, . . . , r". 
Hence the coordinates r* enjoy the following remarkable property: 

Theorem 2 Given a natural n > 2, and two non-negative integers, m G {0, . . . , n — 2} and k, 

k k 

consider the Hamiltonians H^^'^^ = Em,i + X] Cj'^i''^ O'^d H^'^~^ = Em.i + X] CjVi~^\ where Cj 

3=1 ' ' i=i 

are arbitrary constants. 

If k G {!,..., n — m} then the Hamiltonian H^'^^ commutes not only with the 'standard' 

k 

integrals Hmlr^^ = Em,r + J2 CjVr'^\ ^ = 2, . . . , but also with ttj , j = k -\- 1, . . . , n, i.e., are 

i=i 

cyclic variables for -f^m'i^^ /^'^ j = + 1, . . . , n. 

Likewise, if k G {0, . . . ,m} then the Hamiltonian commutes, in addition to Hm,r~^ = 

h 

( '\ ( k ^ 

Em,r + X] Cj^r , T = 2, ... ,71, with 7c j , j = 1, . . . , n — k, i. e., are cyclic variables for ^ ' 
for j = 1,. . . ,n- k. 

Note that, in contrast with the above, the additional integrals of motion for H^f found earlier 
in [20] were quadratic rather than linear in momenta. 
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IV Exact solvability in n dimensions 

Because of the special form ([7]) of the metric Gm in the variables r^,...,r" the existence of 
cyclic variables simplifies solving the Hamilton-Jacobi and the Schrodinger equations for the 
Hamiltonians H^'^^ with m = 0, . . . , n—2 and k = 1, . . . , n—m—l and if^ f ' "* for m = 0, . . . , n—2 
and k = 0, . . . ,m even more than one could expect, especially provided n is sufficiently large. 

Let Vs = —id/dr'^, s = 1, . . . ,n, where i = The following results are readily verified by 

straightforward computation: 

Theorem 3 Fix a non-negative integer m, a natural k, and k constants Cj, j = 1, . . . ,k. Then 
for any natural n such that n > 2k + m the most general common eigenfunction ip of 

/ n—m n \ k 

+ ^ VbP2n-m+l-b I + 'y^.Cj^l^ 

\ a=l b=n—m+l / j=l 

and of Vj, j = k + 1, . . . ,n, with the eigenvalues E and ttj, j = k + 1, . . . ,n, respectively, is 

n 

quasiclassical and reads ip = exp(iS'), where S{r^, . . . , r") = So{r^, . . . , r^) + ^ r%j satisfies 

j=k+i 

the stationary Hamilton-Jacobi equation for H^'i\ and Sq is a general solution of the reduced 
Hamilton-Jacobi equation, a first order linear PDE in k independent variables: 

Errn-n.^i-j^ + Ec,vl^^=e, (9) 
j=i orJ 

where 

_ {n-m-2[{n-m)/2\) ^ 1 ^ 

e — rj — VTn-m+l-jTTj ^ n-m-[(n-m) /2\ ~ ^ 2-^ ^i^2n-m+l-j- 

j=k+l j=n—m+l 

Theorem 4 Fix a natural k, a non-negative integer m > 2k, and k constants Cj, j = 1, . . . ,k. 
Then for any natural n > m the most general common eigenfunction ip of 



^ / n—m n \ k 



. a=l fe=n— m+l / j=l 

and of Vj, j = 1, . . . ,n — k, with the eigenvalues E and ttj, j = 1, . . . ,n — k, respectively, is 

n—k 

quasiclassical and reads ip = exp(iS'), where S{r^, . . . , r") = Sq^t"-"''^^, . . . , r") + ^ r^nj satisfies 

the stationary Hamilton-Jacobi equation for hI^^'~\ and Sq is a general solution of the reduced 
Hamilton-Jacobi equation, a first order linear PDE in k independent variables: 

n f)C k 

E ^2n-m^l-,^+EcM~'^=e, (10) 
j=n-k+l 'Jl j=l 

where 

. _ 1 ^ (m-2[m/2]) , 

^ — E — - TTn-m+l-j'^j — T^jT^2n-m+l-j ^ ^n-[m/2]- 

i=l j=n-[m/2]+l 



Let us briefly outline the integration strategy for and (ITUi) . Consider first. Assume that 
we have found new coordinates z^{r^, . . . ,r''), . . . , z''{r^, . . . ,r^) such that 

This is always possible, but the choice of z^s is not unique and depends on the particular values 
of TTj. Now becomes an ODE in z^ involving z^, . . . , z^~^ as parameters: 

g + Ee,W^'-. (11) 

where V^^ are now considered as functions of 2;'s. 
The general solution of ( fTTI) reads 

^0 = K{z^. • • • , z^-^) - E / c.V^^dz^ + (12) 

i=i 

where is an arbitrary smooth function of its arguments. 

Likewise, assume that we have found new coordinates ^"■-fc+i^^-^-'^+i^ . . . ^ 7-")^ . . . ^ ^n^^n-fc+i^ _ _ _ ^ ^n-j 
such that 



V" d _ d 

jr=n— A;+l 

Then the general solution of fllOl) reads 

= K{z''-^+\ . . . , - E / Cjl^/"-''^ci5" + £5", (13) 

i=i 

where ii' is an arbitrary smooth function of its arguments. 

Thus, the stationary Schrodinger equations for the Hamiltonians H^^'^^ and H^^ f ' •* with 
m = 0, ... ,n — 2 (for -ff^ i ' have an extra condition m > 2k) and arbitrary constants q, 

i = 1, . . . , k, in the space of n dimensions for any n > 2k + m essentially reduce to a linear first- 
order PDEs in k independent variables, and these PDE can be explicitly solved. This is a rather 
surprising result, as only a few potentials for which the Schrodinger equation is exactly solvable 
in the space of arbitrarily high dimension were known so far, cf. the discussion in Introduction. 

It is readily seen that there exists no real value of E for which the eigenfunction 



= exp I i I 5*0 + ^ 

j=k+i 



of Til^'i^ constructed in Theorem [3] with 5*0 given by (fT2|) can have a finite norm 

?/'|W...rfr"~ [ \exp{K{z\...,z''-'))\^dzK..dz''. 
Jq 
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The latter integral obviously diverges for any choice of K because of integration over 2;*^+^,. . . , 
z". Therefore no common eigenfunction ip of 'H^mi^ ^"^^ 'Pjt j = k + 1, . . . ,n, can belong to the 
discrete spectrum of Tif^i- In a similar fashion we can show that no common eigenfunction ip of 
Tiln 1' ^ and Vj, j = 1, . . . ,n — k, can belong to the discrete spectrum of i' \ 



V Example: an anharmonic oscillator 

Consider the Hamiltonian 

k=0 j=k+l 

in the (g, p)-cootdinates. It commutes with 

n 

by construction. 
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A Superintegrability 

The Hamiltonian H is sup erintegr able [20] because it also commutes with the function 



/ = < 



\pI + {q'f, n = 2, 

IpI n = 3, 

Pn, n > 4, 

and hence Fr = {/, Hr}, r = 2, . . . ,n also Poisson commute with H, and all functions in the set 
{Hr, r = 1, . . . ,n, Fs, s = 2, . . . ,n} are functionally independent. 

For n >4 the additional integrals Fr are simply Fr = {pn,Hr} = —dHr/dq^, and we readily 
find that 



i=n~r 

In particular, we obtain 



2 

i=n~r+2 j=2n+2~i—r 



F2=pI/2, (14) 

so F2 is simply a half of square of /, and it is straightforward to verify (see Theorem 5 of [27] for 
a more general result of this kind) that we have 

{pn,Fr} = 0, r = 2, ...,n, {i^,,Fj = 0, r,s = 2,...,n. 
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,l^3 



It can be readily inferred from the above and from Theorem[T]that the quantities {H, I, F^, . . . , Fn} 
are functionally independent and Poisson commute for all n > 2. 

Thus, the Hamiltonian H is maximally superintegrable for all n = 2,3, . . .. For n = 3 we 
have, in addition to H, H2, H3 and J, the following integral K = p^p2 + Q^vil'^ — + {q^Y = 
'^2'^?j — + 3(r^)^/4 which is quadratic in momenta. Note that / and K commute. 

B Quantization and solution of equations of motion 
in the (r, 7r)-coordinates 

In the (r, 7r)-coordinates we have 

1 " 3 1 

j=i 

The quantization is obvious: H goes into the operator 

i=i 

For n > 4 we can look for the common eigenfunctions ip oiTi and Pj, i = 4, . . . , n: 

ni) = E^jj, Vi^p = TTiV, z = 4, . . . , n. (15) 

By Proposition 2, for n > 6 finding requires solving the equation (Q that becomes 

E 7r„+i_fc95o/9r'= - r3 + frV^ - ^{r^ = e, (16) 
fc=i 

where 

. F Y^^ . (n-2K2]) 

k=A 

Eq. (fT6l) is a /irsi order PDE that can be readily solved in full generality. 

Let K{uji,uj2) stand below for an arbitrary (smooth) function of its arguments. 
If 7r„ 7^ then the general solution of f|T6l) reads 



+K 



7r„ vr„ 
If 7r„ = but 7r„_i 7^ then the general solution of ( fT6l) is 

(^1)3^2 3,.1(,.2^2 ^^_^(^^2^2 (^3_^^)^2 / (r^7r„_i - 7r„_2r^) 

= -7: z TT^ \ h A r , 

2vr„-i 47r„_i ^Tr^^;^ 7r„_i \ 7r„_i 
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Finally, if 7r„ = 7r„_i = but 7r„_2 7^ then the general solution of (1161) has the form 

^ (r^)^ — Sr^r'^r^ + (r^)^r^ + 2er^ _|_ i 2-j 
2vr„_2 

By Proposition 2 for n > 6 the most general common eigenfunction ip oi Ti and of Vi, i 
4, . . . , n is 

n 



t/- = exp [ i [ 5*0 + ^ 

i=4 



where 5*0 is given above. This eigenfunction is not square integrable for any choice of K and for 
any values of tTj. 

An integral of the stationary Hamilton-Jacobi equation for H is given by 

n 

S = So + y^TTj-r-^. 

i=4 

C Special cases: n < 5 

We list below the exact solutions for the corresponding Schrodinger and Hamilton-Jacobi equa- 
tions. None of the eigenfunctions listed below is square integrable. 

Case 1: n = 2 

Here = and there is a pair of commuting operators H and I. In the (r, 7r)-coordinates 
we have 

n = V1V2 + 2rV2 - (r^)V2, J = Vi/2 + {r^f. 

and their common eigenfunction ijj satisfying Tiip = Eijj, Iijj = Xiijj has, up to multiplication by 
an arbitrary constant, the form 

= r + vW^) exp (v/2((ri)2-A0(r^ - (r'f/U - X,/6)) 

We actually have two such eigenfunctions, as a/ (r^)^ — Ai is a two-valued function. 

The eigenfunction ip is not quasiclassical: we have a complete integral of the stationary 
Hamilton-Jacobi equation for H of the form 

S = -i\n^P + ^{{r' f - Ai)^/^(-2Ai - {r')' + 12r') - ln{ir'f - Ai), 



or equivalent ly. 



^^^arctan^ ^ 



^ \{{r^)2 - J 12 

11 



i-(12r2-(rY-2A0(2(Ai-(ri)2))^/2. 



Case 2: n — 3 

Now we have a triplet of commuting operators H, I — 7:^/2 — r^, and {F2 — Vl/2 — ?,r^Vz/2 — V2 
and their common eigenfunction such that 'H%1) — Eip, lip — Xiip and ^2'0 — ^2^^ reads 

^ = (^1 ^\y/2 exp(iV2(Ai + r'y/'r' + i(-A2 + V2{r' + A^V^A^r^ 

-i(-A2 + ^(r^ + AO^/^AOAir^ + IT^^' + A0^/^(-10A? + 14£; - 7A^)) 

If we take another triple of commuting operators Ti,! — 7r|/2— r^, and /C = 7^27^3— r'^+3{r^y /A 
then their common eigenfunction such that Hip — Eip, lip — Xiip and fCip — vip reads 

^ - FT^^^^(^^^'^ ^^^^'^^^ - ^/2(A. + .T/- ("^^^^^^ 

+ \r\r'f + ^(-7^ + 2AD(r^)^ - \{UE + 2A? - 7A,^)r^ - (r^)V2 - ^r^ 



-(-2A2i/+|At + 2AiE+ii/2)j^. 



In the latter case the eigenfunction ip is almost quasiclassical: we have a complete integral of 
the Hamilton-Jacobi equation for H of the form 

,S=-iln^ + ^ln(r^ + A^). 

In both cases we again actually have two eigenfunctions as {r^ + Ai)^/^ is two- valued. 

Ceise 3: n = 4 

The common eigenfunction of 7^4, Ti. and i = 2, 3, with the respective eigenvalues 774, A2,A3 
reads 



= exp(i7r4r^ + ^(27rl + 2 + 37rl)(r^)^ + ^A^Itt^ + 2)(r^)3 - ^(37r| + 2)(r^)V2 



(-3A^ + 2^lX,){T'f - %rVi + ^rh' + -i.(-2A27r42A3 + A^ + 27r4^£;)ri 



47rr ■ ^ TT^ • ■ 7r4' • ' 271 



27r4 27r| 

This eigenfunction is quasiclassical: 5* = — iln-^ is a complete integral for the stationary 
Hamilton-Jacobi equation for H. 
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Case 4: n — 5 

The common eigenfunction of Vi, i — 4,5, H and Ti, i — 3,4, with the respective eigenvalues 
7r4,7r5, E, A3,A4 reads, up to an overall constant factor, 

--^(-S^TT^ - + 37r^ + 87r4A47r5 - 4'Kl\^)r^ + -^(r^)^ 

oTTg Z7r5 

+ 2^(2A47r5 - 27r4A3 + Trfjr^ - ^(tt^ - 2A3)r3) 

Unlike the previous case, this wave function is not quasiclassical: there is a complete integral 
of the stationary Hamilton-Jacobi equation for H of the form 
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